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I .  INTRODUCTION 


The  design  of  radar- stealthy  platforms  has  become  an 
important  engineering  problem  over  the  last  several  decades. 
Although  the  principles  of  radar  stealth  are  well  known,  radar 
cross  section  (RCS)  prediction  and  reduction  are  still 
difficult  tasks.  The  development  of  increasingly  sophisticated 
detection  systems  threatens  to  reduce  the  mission 
effectiveness  of  many  types  of  weapons  platforms.  Much  effort 
is  now  spent  on  reducing  radar  detectability  by  reducing  the 
RCS  of  the  platform. 

Four  basic  techniques  are  employed  in  radar  cross  section 
reduction  (RCSR) : 

1.  shaping, 

2.  surface  material  selection, 

3.  active  cancellation,  and 

4.  passive  cancellation. 

Each  method  has  advantages  and  disadvantages  and  the  final 
selection  is  a  compromise.  Shaping  is  one  of  the  most 
important  ways  of  controlling  the  radar  echo,  but  usually  RCSR 
at  one  viewing  angle  is  accompanied  by  an  enhancement  at 
others.  Radar  absorbing  materials  (RAM)  can  also  be  used,  and 
the  RCS  reduction  for  this  case  is  achieved  by  the  dissipation 
of  energy  inside  of  the  material.  The  platform  performance 
(e.g.,  aerodynamic  efficiency)  generally  decreases  due  to 
added  weight.  Surface  maintenance  requirements  are  also  more 
demanding  than  those  for  conventional  materials.  Active  and 
passive  cancellation  are  limited  methods  and  can  only  be 
applied  to  narrow  frequency  bands  and  spatial  regions. 

When  choosing  materials  it  is  necessary  to  establish  a 
relationship  between  the  electrical  characteristics  of  the 
material  and  the  target's  scattered  field.  The  electrical 
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properties  of  a  material  are  usually  expressed  in  terms  of  the 
complex  permittivity  and  permeability  (e  and  n  respectively) . 
Finding  the  scattered  field  when  e,  fi  and  the  shape  of  the 
target  are  known,  is  straightforward,  although  not  always 
easy.  RCS  synthesis,  the  inverse  problem  (i.e.,  finding  £  and 
/x  for  a  given  RCS  pattern  and  body  shape)  is  more  difficult 
because,  in  general,  a  complicated  set  of  coupled  integral 
equations  must  be  solved.  Because  of  the  complexity  of  the 
synthesis  procedure,  RCS  reduction  using  RAM  has  been  based  on 
intuition  and  experience. 

An  RCS  synthesis  method  for  an  impedance  surface  can  be 
based  on  the  solution  of  the  scattering  integral  using  a 
method  of  moments  (MM)  technique.  The  (unknown)  surface 
impedance  is  expanded  in  a  series  of  basis  functions  with 
unknown  coefficients.  In  the  most  general  synthesis  case,  both 
electric  and  magnetic  currents  must  be  considered.  When  the 
surface  is  planar,  an  expression  for  the  currents  can  be 
obtained  by  applying  an  approximate  set  of  boundary  conditions 
at  the  surface  [Ref.  l] .  In  general,  a  coupled  set  of  integral 
equations  must  be  solved  [Ref.  2].  The  solution  yields  the 
surface  impedance  at  every  point  on  the  surface  of  the  body 
for  a  specified  scattered  field  (or  RCS)  pattern.  Once  the 
surface  impedance  is  known,  the  electrical  properties  of  the 
material  can  be  determined  under  certain  conditions  via  the 
refractive  index  of  the  material  [Ref.  3]. 

A  resistive  sheet  is  a  special  case  of  an  impedance 
surface  since  it  does  not  support  magnetic  currents.  It  is  of 
great  importance  for  RCS  reduction  purposes  because  it  is  used 
to  control  travelling  wave  effects.  In  this  thesis,  the 
synthesis  procedure  described  above  is  applied  to  a  resistive 
sheet,  in  which  case  the  integral  equations  simplify 
considerably. 

Chapter  II  is  a  review  of  the  concepts  and  formulas 
required  in  the  formulation  and  solution  of  the  problem.  The 
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reflection  and  transmission  coefficients  for  a  thin  resistive 
surface  are  also  derived  from  the  boundary  conditions.  In 
Chapter  III  the  appropriate  equations  are  derived  and 
expressed  in  matrix  form  for  the  monostatic  and  bistatic 
cases.  Chapter  IV  deals  with  the  computer  implementation  of 
these  equations  and  the  evaluation  of  the  obtained  data. 
Finally,  Chapter  V  concludes  with  a  discussion  on  the 
benefits,  concerns,  and  recommendations  for  future  work. 
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II.  THEORETICAL  BACKGROUND 


In  this  chapter,  a  few  concepts  and  equations  used 
throughout  the  thesis  are  introduced. 


A.  RESISTIVE,  CONDUCTIVE  AND  IMPEDANCE  SURFACES 


Thin  layers  of  lossy  materials  are  of  special  interest  in 
RCSR  design  as  they  have  lower  RCS  than  the  corresponding 
perfect  electric  conductors  [Ref.  3] .  A  mathematical  model  of 
such  a  layer  is  an  electrically  resistive  sheet.  It  can  be 
visualized  as  a  thin  sheet  of  highly  conducting  material  whose 
permeability  is  that  of  the  surrounding  medium  (free  space) . 
Theoretically,  it  is  composed  of  an  electric  surface  current 
of  infinitesimal  thickness  whose  total  strength  is 
proportional  to  the  tangential  electric  field  at  its  surface. 
The  sheet  is  characterized  by  a  jump  discontinuity  in  the 
tangential  components  of  the  magnetic  field  across  the 
surface,  but  no  discontinuity  in  the  tangential  electric 
field.  Therefore,  a  resistive  sheet  does  not  support  magnetic 
currents.  Its  properties  are  completely  specified  by  its 
resistivity,  Rs  in  ohms  per  square.  This  unit  is  derived  from 
the  basic  definition  of  the  resistance  of  a  rectangular  cell 
with  length  1,  and  cross  section,  S.  For  the  special  case  of 
length  equal  to  width  (l=w) 


=  — 4—  =  -^r  [-Si 

s  oS  o  tw  at  □ 


(1) 


where,  a  is  the  conductivity  of  the  material  in 
(Siemens /meter)  and  t  is  the  thickness  of  the  cell. 

The  boundary  conditions  at  the  two  sides  of  the  sheet 
follow  from  the  usual  boundary  conditions  of  electromagnetics 
[Ref.  2]  .  For  the  limiting  case  of  t  approaching  zero,  a 
increases  in  such  a  manner  that  Rs  is  finite  in  the  limit. 
However,  when  Rs  =  0  the  sheet  is  perfectly  conducting,  and 
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when  Rs-°°  it  is  perfectly  transparent.  Thus,  the  boundary 
conditions  for  the  resistive  sheet  are  [Ref.  4] 

A//Un  =  Js  =  otTun  =  ^ 

Ks 

Etan(  +  )  =*tan(")  =ESJS 

where  the  signs  refer  to  the  upper  (positive)  and  lower 
(negative)  faces  of  the  sheet,  and  Js  is  the  total  surface 
electric  current  supported  by  the  surface. 

A  general  vector  form  of  the  boundary  conditions  are 
[Ref.  1] 


fixE(  +  )  -fixE{~)  =  0 

(2) 

A  xH{  +  )  -AxH(-)  =JB 

(3) 

fix{fixE{±)}  =  ~RSJS 

(4) 

where  h  is  the  unit  vector  normal  to  the  sheet  and  directed 
towards  the  positive  (plus)  side  of  the  sheet  and  Js  is  the 
total  electric  current  supported. 

The  electromagnetic  dual  of  the  resistive  sheet  is  the 
magnetically  conductive  sheet  having  a  conductivity  Rm  in 
S/sq.  Such  a  sheet  has  a  jump  discontinuity  in  the  tangential 
electric  field  but  none  in  the  magnetic  field,  and  supports 
only  magnetic  currents.  It  simulates  a  thin  layer  whose 
permeability  differs  from  that  of  the  surrounding  medium  (it 
would  be  difficult  to  realize  this  type  of  sheet) .  When  Rm= 0 

the  sheet  acts  like  a  "perfect  ferrite"  with  infinite 
permeability  and  when  Rm  =  eo  it  no  longer  exists.  The  general 
boundary  conditions  for  such  a  sheet  are 
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with 


nxH (  +)  -M(-)  =0 
J3x£(+)  -fixE(-)  =-jm 


fix{fixH{±)}=-RmJm 

where  Jm  is  the  magnetic  current  supported  and,  the  convention 

followed  for  the  signs  is  the  same  as  in  the  resistive  sheet 
case . 

Note  that  an  impedance  surface  is  a  simulation  of  a  thin 
layer  whose  permittivity  and  permeability  both  differ  from  the 
surrounding  medium.  Therefore,  a  combined  sheet  consisting  of 
resistive  and  conductive  ones  is  equivalent  to  an  impedance 
surface.  Although  the  two  sheets  are,  in  general,  coupled, 
with  each  affecting  the  scattering  from  the  other,  decoupling 
occurs  when  the  sheets  lie  in  a  plane.  It  is  therefore 
sufficient  to  solve  the  problem  for  the  simpler  case  of  the 
resistive  sheet  and  use  the  duality  principle  to  get  the 
solution  for  the  other.  The  solution  of  the  combination  sheet 
can  follow  by  adding  the  two  solutions  [Ref.  1] . 

B.  REFLECTION  AND  TRANSMISSION 

The  reflection  and  transmission  coefficients  for  a 
resistive  film  can  be  derived  from  the  boundary  conditions  in 
Equations  (2)  ,  (3)  and  (4)  and  the  geometry  shown  in  Figure  1. 
The  subscripts  i,  s  and  t  refer  to  incident,  scattered  (or 
reflected)  and  transmitted,  respectively.  The  two  principal 
polarizations  (parallel  and  perpendicular)  are  investigated 
separately.  Specular  reflection,  and  transmission  at  the  angle 
of  incidence  is  assumed,  because  the  sheet  is  infinitely  thin. 


7 


Figure  1:  Reflection  and  transmission  from 
a  resistive  sheet. 


1.  Parallel  Polarization 

For  parallel  polarization,  the  electric  field  of  the 
incident  wave  lies  in  the  plane  of  incidence 

5i=0£’eie'iPi'i’.  (5) 

Without  loss  of  generality,  it  can  be  assumed  that  the  plane 
of  incidence  is  the  x-z  plane.  Figure  2  shows  a  cut  through 
the  (J)  =  0  plane. 


Figure  2:  A  cut  of  Figure  1  in  the  plane  of 
incidence . 
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where  rp  is  the  reflection  coefficient  for  parallel 
polarization.  Using  the  notation  introduced  in  Section  A  and 
applying  the  boundary  conditions  on  the  electric  field 
intensity  at  z=0 
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E(+)  =  (E^ES)  |z  =  0, 

(12) 

E(~)  =  ^U- 

Similarly  for  the  boundary  conditions  on  the  magnetic  field 
intensity 

H(  +  )  =  {Hi+Hs)  \Z=Q  , 

(13) 

H(~)  =  |z.0. 

Inserting  Equations  (6),  (8),  (10),  and  (12)  into  (2)  gives 

Tp  =  l+Tp.  (14) 

Using  Equations  (7),  (9),  (11)  and  (13)  into  (3)  gives 

(i-r  -t  )  e'J'p[“x£in(^,]  .  (15) 


Using  Equation  (4)  and  Equation  (15! 


TpCOS  ( 0i )  (l-rp-Tp)  . 


Solving  Equations  (14)  and  (16)  for  the  reflection  coefficient 
yields 

_  -T^coste,) 

P  2i?s+riC)Cos(ei)  1  * 

and  for  the  transmission  coefficient 


2  Rs 

P  2J?£  +  TloCOS(ei) 


10 


2 .  Perpendicular  (Normal)  Polarization 

For  perpendicular  polarization  a  similar  procedure  is 
followed  for  deriving  the  reflection  and  transmission 
coefficients.  In  general  the  incident  field  is  of  the  form 


E1=$E0e 


(19) 


Again,  without  loss  of  generality,  the  x-z  plane  is  considered 
as  the  plane  of  incidence.  Expressions  for  the  incident  waves 
are 

Ei=-fiE  [-^sinfBi)  -zcos(6i)]  ,  (20) 


[-Scosfe,)  *i?sin(e2n  .  (21) 

Vo 

Denoting  T„  and  r„  the  transmission  and  reflection  coefficients 
respectively  for  normal  polarization,  the  scattered  and 
transmitted  fields  are 

Et  =  -9T  E  e"lP  [-xsint0*)  -zcos^)]  ,  (22) 


H ' 


Tn—  [-xcos  (0i) 

Tlo 


Ee=-?rnE0e 


-jp  [-jeeinOi) +zcos  (Bi)  ] 


(24) 


and 
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(25) 


Hs  =T  E° 


n  Ucos  (Oj)  +2 sin  (0i)  ]  e 

'l  o 


iP  [-xsin  (0j)  t-zcos  (0^)  ] 


Applying  Equations  (20),  (22),  (24)  and  (12)  to  (2)  gives 

Tn  =  l^Tn.  (26) 

To  obtain  the  surface  current,  Equations  (21),  (23),  (25)  and 

(13)  are  inserted  into  (2)  yielding 

Js  =  9—c 05(0,)  (rfl  +  ra-l)  e-j*l-XBin{e<)]  .  (27) 

*1  o 


Using  Equation  (27)  in  the  boundary  condition  (Equation  (3)) 
gives 


~  COS  (0j)  (l-T„-rj  . 

I  o 


(28) 


Solving  Equations  (26)  and  (28)  for  the  reflection  and 
transmission  coefficients  gives 


*1  o 


2Rscos  {Qd)  +T|0 


(29) 


and 


2i?Bcos  (0i) 
2Rscos  {dd)  +r)0  ‘ 


(30) 


C.  RADIATION  AND  SCATTERING  EQUATIONS 

In  general,  electric  and  magnetic  currents  are  the 
sources  of  radiated  or  scattered  fields.  For  the  case  of  a 
resistive  sheet  the  magnetic  currents  are  zero.  The  formulas 
simplify  further  when  only  the  far  zone  fields  are  of 
interest . 
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The  formulas  for  the  far  -  field  calculation  are  derived 
by  Balanis  [Ref.  5]  and  are  given  in  this  section  for  the 
reader's  convenience.  They  are  specified  for  the  case  of 
currents  in  the  x-y  plane  of  a  Cartesian  coordinate  system  as 
shown  in  Figure  3 . 

The  spherical  components  of  the  electric  field  intensity 

are 

Er“  0  , 


*e  = 


J  P  e~J 

4  r 


Tlo-W'e/ 


(31) 


j  p  e  jPr 


4  7t  X 


where 


(32) 


Ne=JJ  [  Jxcos(0)  cos  (<!>)  +Jycos  (6)  sin  (<|>)  ]  e^^^dx' dy'  (33) 


N^ff  [-<Jxsin((J))  +Jycos(<|))]  e^r,C0B* dx'dy' ,  (34) 


and 

r  ^osilf  ^x'sin  (0)  cos  (4>)  +y/sin  (0)  sin  (<j>)  =x' u+y'v= g.  (35) 

For  the  observation  point  appearing  in  Equation  (35)  the 
direction  cosines  are  defined  by 

u  =  sin  (0)  cos  (<j>)  , 

v=sin(0)  sin(<J>)  ,  (36) 

w=cos  (0)  . 

In  Chapter  III  the  currents  derived  in  Section  B  will  be 
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used  in  the  radiation  integrals  (Equations  (31)  and  (32))  to 
obtain  the  scattered  fields.  The  currents  are  functions  of  the 
surface  resistivity  through  the  reflection  (or  transmission) 
coefficients.  Thus,  for  the  synthesis  problem  Ee  and  in 
Equations  (31)  and  (32)  are  known  but,  the  currents  Jx  and  Jy 
(functions  of  Rs)  are  unknown. 


Figure  3:  Coordinate  system  for  scattering 
from  a  planar  rectangular  plate. 

All  the  formulas  assume  that  the  current  is  not  affected  by 
the  presence  of  the  edge.  This  turns  out  to  be  a  valid 
assumption  for  resistive  surfaces  of  finite  size  because  the 
edges  are  "softer"  than  perfectly  conducting  edges. 
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III.  FORMULATION  AND  SOLUTION 


The  scattering  from  a  planar  rectangular  resistive  sheet 
with  surface  resistivity  Rs  is  investigated  for  the  monostatic 
and  bistatic  cases.  Using  Equation  (4),  an  estimate  for  the 
current  distribution  is  obtained,  and  the  total  scattered 
field  is  calculated.  To  verify  the  synthesis  equations,  the 
calculated  field  is  used  to  initiate  the  synthesis  procedure 
for  the  surface  resistance. 

The  geometry  under  consideration  is  still  that  shown  in 
Figure  1.  The  sheet  lies  in  the  x-y  plane  and  is  illuminated 
by  a  plane  wave  whose  electric  field  contains  both  parallel 
and  perpendicular  components 

E1  =  {6Eq+$E^)  <37> 


where 

-j'Pyf  =  j  P  [x  sin(0i)  cos  (<1^)  +y  sin(0i)  sin  (4^)  +z  cos  (0J  ] 

=  j'P  [x  Uj+y  v±+z  w£]  . 

On  the  surface  of  the  sheet  [z  =  0  and  r  =  r') 

-j$i-r=-j$i-r,=j${x,u1+y'vi)  =jfih.  (38) 

The  formulation  as  well  as  the  analysis  of  the  problem  is 
different  for  the  bistatic  and  monostatic  cases  and, 
therefore,  each  is  considered  separately. 

A.  BISTATIC  CASE 

Equation  (4)  with  h  =  2  becomes 

2 x  [2xE{-))  =  -RSJS  (39) 


where 
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(40) 


E(-)  =  EC  =  (GTpEj+$TnEj)  e^h  . 

The  cross  products  are 

2x&  =  -2  c os  (0)  sin  (<J>)  +j>cos  (0)  cos  (<{>) 
2 x$  =  -£cos  (<}>)  -^>sin(<})) 


and  the  triple  products 

2 x  (2 x0)  =  -  (2cos  (0)  cos  (<J>)  +ycos  (0)  sin  (<J>)  ) 

and 


2 x  (fx$)  =  jfsin  (<J>)  -ycos  (<J>)  . 


Then 


2x  {2xE{~)  )  -{-2 [cos  (Oj)  cos  3^ -sin ((fjJ  Tn e£ ] 

~9  [cos  (0j)  sin  (<|)i)  rpJE6+cos  (4^)  Tn£^]}eJ'ph. 

Substituting  in  the  explicit  expressions  for  the  transmission 
coefficients  (Equations  (18)  and  (30)),  into  (40),  Equation 
(39)  gives  the  components  of  the  current  on  the  surface  of  the 
sheet 


f  COS  (Oj)  COS  (ifrj)  Eq1 
{  2i?s+ri0cos  (0i) 


sin  ( 4>j )  cos  (0,)  E$ 

2  Rs COS  (0j)  +t|  D)  , 


e^h , 


(41) 


-  2 


cos  (0j)  sin  (4^)  Eq 
,  2-Rs+t)0cos  (0J 


+ 


COS  (fa)  COS  (Oj)  Eq 
2  Rb  COS  ( 0_j )  +T)0  ; 


eip/7  _ 


(42) 


Now  that  the  currents  are  known,  the  scattered  fields  can  be 
calculated.  Because  of  the  complexity  of  the  formulas, 
parallel  and  perpendicular  polarizations  are  treated 
separately.  Besides,  any  arbitrary  polarization  can  be 
decomposed  into  a  sum  of  these  two  orthogonal  components . 
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1.  Parallel  Polarization 

When  the  field  is  parallel  polarized  =0  and  Equations 
(41)  and  (42)  simplify  to 

J^2£8VP'’cos(9i)cos(4)i)  2J^os(6j)  .  <«> 


J  =2Eq2  e^hcos  (0i)  sin(4>i) 


2J2s+TloCOS(0i) 


Using  equations  (43)  and  (44)  in  (33)  gives 

Nq  =  2  Eq  cos  (0)  cos  (0_f)  cos  ((J)-^)  Ix  (45) 

where  a  new  quantity  has  been  defined  for  convenience 


r  r  (h+g) 

Ji-ff - - - r5— -dx'  dy' 

J  Js  2Rs+T\0COS  (0j) 


The  resistivity  is  a  function  of  position  in  general  and, 
therefore,  it  cannot  be  taken  out  of  the  integral.  In  order  to 
reduce  the  scattering  equations  to  closed  form  expressions,  it 
is  assumed  that  the  resistivity  profile  can  be  accurately 
represented  by  a  step  approximation.  The  resistivity  function 
is  represented  by  a  series  of  basis  functions  (pulses)  with 
unknown  coefficients 


Rgix'.y')  =^2^2rm,fmn(x,ly/) 


where  M  and  N  are  the  number  of  patches  in  the  y  and  x 
directions  respectively.  The  unknown  expansion  coefficients 
are  { rmn }  and  the  expansion  functions  are  two-dimensional  steps 
or  "pedestals" 
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(48) 


x'<  y'  e  smn , 

else . 

The  geometry  is  defined  in  Figure  4.  The  center  of  the  mnth 
patch  is  at  (xm„,  ymn)  and  each  patch  has  dimensions  Ax  by  Ay. 


fmn  y') 


i  , 

o  , 


Ay 


Figure  4:  Division  of  the  sheet  into 
rectangular  subdomains. 


Applying  the  series  representation  (47)  to  (46)  gives 

rr  eiP  (h  +  g) 

1ocos(ei)dxdy 

r  r  eJ  P 

J  J  s  "  " 


dx'dy' 


w=l  n= 1 


M  N 


2EErmn^+Tioc°s(ei) 

eJp  {h+g) 


For  the  miih  term  of  the  sum  the  integral  is 


rr  eJ'P  .  . 

//  -z - T5-rdx/dy/  = 

J  2rm  +n,cos  (0, 


mn 


i-T]0COS  (0j) 

— - - - t q-t-  f  f  eJ'k(h+g)  dx'dy ' . 

2rmn  +  TloCOS  (6i)  •'•'San 


(49) 
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Evaluating  the  integral  in  (49)  yields  a  closed  form 
expression  in  terms  of  sine  functions.  Thus  Ix  becomes 


Ii  =  £  amne: 


m=l  n=l 


where 


_  sin  [p  (u+Uj)  Ax/2 ]  .  r  p  (u+Uj)  Ax] 

5"  =  fT^^TA^72  ~sinc[  2  J  ' 


sin [p  (v+vA  Ay/ 2]  .  P  {v+vj  Ay 

S*-  IHv+v/a y/2—-SlnCl - 2 - ]  ' 


2*mn+1\ocosWi)  ' 


Using  Equation  (50)  allows  Equation  (45)  to  be  written  as 

M  N 

Nq  =  2Eq  cos  (Gj)  cos  (0)  cos  (<| AxAySuSv^2  £  amne3^Pim. 


Finally,  the  complete  expression  of  the  6- component  of 
the  scattered  electric  field  for  a  parallel  polarized  incident 


field  is 


where 


Eqq  =  CdEqC os(0i)  cos  (0)  cos  (4)-^) 
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(56) 


-j'Ptlc,e_J'pr 

°  2  7t  r 

The  first  subscript  on  Epqs  denotes  the  polarization  of  the 
scattered  field,  whereas  the  second  subscript  denotes  the 
polarization  of  the  incident  wave. 

By  the  definition  of  the  radar  cross  section  for  a  three- 
dimensional  target 


o=l im  [4t ir2  1-  ,1  ] 
r-  \E2\2 


one  obtains  for  this  case 


°ee  “ 


(rioP: 


[cos  (0i)  cos  (0)  cos  ( <j> ~4>j: )  Tj 


To  calculate  the  0- component  of  the  scattered  electric 
field,  E^,  Equation  (32)  is  used  along  with  Equations  (34)  , 
(43)  and  (44).  Following  the  same  procedure 


£"4,0=  -C0EqCOS  (0i)  sin  (<J)-(j)i)  I1 


and  the  corresponding  radar  cross  section  is 


°4>0  ~ 


(iloPi 


[cos  (©j)  sin  ((J)-^)  Jj  2  . 


2.  Perpendicular  Polarization 

When  the  incident  field  is  normally  polarized  (Ej  =0)  , 
the  current  components  in  Equations  (41) ,  (42)  reduce  to 

^-2^e^sln(W  costf,)  2Jg>eos1(flj)  ,  (61) 
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(62) 


Jy  =  2E^e^hcos  ( <j>i )  cos  (0i) 


2 i?scos  (Oj)  +rj ( 


Using  Equations  (61)  and  (62)  in  (33)  gives 


Nq-2 cos  ( 0 )  cos(0i)  sin  ((j)-^)  I2 


where 


rr  ejk(h+g) 

I2=ff  O  E>  rinic  I  ft  TTZ-dx'dy'  . 

J  J  s  2  RsCOS  ( 0  ^ ) 

Following  the  procedure  used  for  parallel  polarization  one 
obtains 


J2  =  AxA  ySuSv  EE  mn  ® 


with 


^  2  TjnnCOS  (0j)  +T10' 


Finally,  using  Equation  (63)  in  (31)  gives  the  6 -polarized 
component  of  the  scattered  field 

EqI  =  C0E$  cos  (0)  cos  (0j)  sin(<|)-<j)i)  J2  .  (66) 

The  corresponding  radar  cross  section  is  then 


(TloP)2 


[cos  (0)  cos  (0i)  sin  (4>~4>i)  I2] 2  ■ 


To  determine  the  0- component  of  the  scattered  field,  Equation 
(34)  is  used,  yielding 


•Afy  =  2  Eq  cos  (0i)  cos  (<|>-<i>i)  J2 


and  then 


E^  =  C0E^cos  (0j)  cos  (4)-^)  J2  . 
The  corresponding  radar  cross  section  is 
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(69) 


(f)oP) 


[cos  (0i)  cos  ((J)-^)  J2] 2  ■ 


3 .  Arbitrary  Polarization 

In  general,  the  incident  wave  is  composed  of  both 
parallel  and  perpendicular  components.  Thus  it  is  necessary  to 
combine  the  results  of  Subsections  1  and  2.  The  total  8  and  $ 
components  of  the  scattered  field  intensity  are  given  by 

Eq  =  Eq%  +  Ee£^ 

(70) 

=  C0cos  (0)  cosiQ^  [EqCos  )  J1+£'^sin(<j)-(j)i)  J2] 


Ef  -  jEufo  +  Emm 


=  -Cocos(0i)  [Eq1  sin  I1-E^1cos  (<J)-<J)i)  J2)  . 

At  this  point  it  is  possible  to  cast  Equations  (70)  and 
(71)  in  matrix  form 

_  ^ee  si  j72) 

S-, 


where,  using  Equations  (50)  and  (64) 

AqG=Eq  C0AxAycos  (0i)  cos  (0)  cos  ((J)-^)  SUSV 


Aqq-Eq  C0  Ax  Ay  cos  (0J  cos  (0)  sin  (<J)-<l)i)  SUSV 
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At>e  =  ~Ee  CaAxAycos  (0j)  sin  ((j)-^)  SUSV 


(75) 


AH  =  C0AxAycos  {QJ  cos  (0-4^)  5U5V.  (76) 

The  terms  S2  and  S2  in  the  right  hand  side  of  Equation  (72) 
represent  the  double  summation  terms  in  (50)  and  (64)  for  each 
scattered  field  point,  namely 

SrEEv,,P-  (77) 


The  resistivity  is  embedded  in  the  terms  Sj  and  S2. 


B .  MONOSTATIC  CASE 

For  the  monostatic  case  (0=0;  and  6  =  6 ;)  and  Equations  (70) 


and  (71)  reduce  to 


Eq  =  CqEq  cos2  (0)  Ix 


Ef  =  CnEt cos  (0)  J2 


respectively.  Equation  (53)  can  be  written  now  as 


2r^+ti0cos  (0) 


and  Equation  (65)  as 


(82) 


Now  it  is  possible  to  set  up  a  linear  system  of 
equations.  In  fact,  as  a  6 -polarized  incident  electric  field 
creates  only  0-polarized  scattered  electric  field,  and 
similarly  for  a  ^-polarized  incident  field,  an  equation  of  the 
form  of  (72)  contains  only  one  component  of  the  scattered 
field.  These  equations  were  implemented  in  MATLAB  and  the  code 
listings  are  given  in  the  Appendix. 
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C.  SYNTHESIS  FORM  OF  THE  SCATTERING  EQUATIONS 

For  RCS  synthesis,  the  scattered  field  components  derived 
in  the  two  previous  sections  are  known  quantities  and  the 
expansion  coefficients  { amn }  and  must  be  determined. 

Expressing  the  scattered  field  equations  in  appropriate  matrix 
form,  it  is  possible  to  solve  for  the  expansion  coefficients 
using  standard  matrix  techniques.  To  achieve  this  it  is 
convenient  to  use  a  single  index  to  count  the  patches  rather 
than  the  double  summation.  Let 

M  N  MxN 

(87) 

m=l  n= 1  k- 1 

where  cmn  can  be  either  amn  or  bmn. 

There  are  MxN  unknowns  and  therefore  MxN  equations  are 
needed  to  solve  for  {amn}  or  {bmn} .  Because  of  the  periodicity 
of  the  sine  functions  (Equations  (51)  and  (52))  in  the 
direction  cosines  u  and  v,  it  is  advantageous  to  work  in 
direction  cosine  space.  As  shown  in  Figure  5,  the  region 
defined  by  -Is  u,  v  <1  is  divided  into  the  same  number  of 
subsections  as  there  are  patches  on  the  sheet.  The  centers  of 
these  subsections  (uk,vk)  form  a  grid  of  points  at  which  the 
scattered  field  components  are  specified,  yielding  MxN  values 
for  each  field  component.  The  area  in  the  unit  circle  is  the 
visible  region  while  exterior  points  have  no  physical 
significance . 
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1 


V 


(uk  ,Vk) 


Figure  5:  Division  of  direction  cosine 
space  into  subsections. 


The  values  of  Ef  and  are  used  to  form  a  column  vector 
of  length  MxN.  Furthermore,  at  each  of  the  patches,  cos(Ok), 
sin  (4>k-<i>i)  ,  cos  (4>k-<t>j)  can  be  expressed  in  terms  of  (uk,vk) 
according  to  the  following  relationships 

ujt  +  vj^sin^e*)  , 

—  =  tan  ($*)  . 
uk 

Inverting  these  formulas  gives 

cos  (6*)  =^1-  (ua2  +  v|)  , 

(88) 

=  arctan  ( — )  . 

uk 

In  the  direction  cosine  space  the  terms  cos (8),  sin ($-&)  , 
cos  (cp-fy)  as  well  as  the  product  SBSV  can  be  written  in  the  form 
a  column  vector  of  length  MxN.  Denoting  column  vectors  with  an 
underline 
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(89) 


pyfc  =  cos  (Qk)  , 
[Z2]k  =  cos(4)Jc-<l)i)  , 
[Z3]jc  =  sin(<l)Jt-<l)i)  , 

FJ*=5«(U*)  Sv[v*] 


We  can  further  define  the  column  vectors  X^  and  X^  as 


MxN 

E 

V=1 


a  giP  Pv(Ujc>vk) 


(90) 


and 


MxN 


-E 


bve 


J  P  Pv  ( uk,  vk) 


V  —1 


(91) 


The  vectors  and  Xf,  can  also  be  written  as 


Xa=Ba  (92) 

and 

Xb  =  Bh  (93) 

where  a  and  b  are  column  vectors  of  length  MxN  and  follow 
from  equations  (53)  and  (65)  respectively.  The  matrix  B  has 
dimensions  (MxN)  by  (MxN)  and  is  defined  as 

=  ei  P [xv  ( +yv  ( vk+v±)  ]  (94) 


where  k,  v  =  1 ,  .  .  .  ,  MxN. 

With  the  above  matrices  defined,  Equation  (70)  can  be 
written  as 
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(95) 


E*  =  qc os  (0,)  q  *X4(E6iX2*Xa+EjX2*Xb) 

where  *  denotes  an  elementwise  product  and  C1  =  C0AxAy. 
Equation  (71)  becomes 

M%  =  -qcos  (0,)  X4  *{Ee1X3  *Xa-EjX2  *Xb)  .  (96) 

For  the  monostatic  case  Xa  =  Xb  =  X  and  Equation  (83)  becomes 

Eq  =  C± Eq  X  (97) 

and  Equation  (84) 

E^C.Ejx^X.  (98) 

For  this  case  u,=u  and  v,.=  v,  and  the  matrix  B  is  defined  now  as 

[B]kv  =  e2j^Uk^v^  (99) 

In  all  of  these  equations  the  goal  is  to  solve  for  the  vectors 
So.  &  or  X/  which  can  easily  be  done  using  elementwise 
division.  Once  these  vectors  are  known,  it  is  possible  to 
solve  for  a  or  b  of  Equations  (92)  and  (93) . 
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IV.  COMPUTER  IMPLEMENTATION  AND  DATA  ANALYSIS 
A.  GENERAL  APPROACH 

In  this  chapter  the  surface  resistivity  synthesis 
procedure  is  simulated  and  computational  accuracy  issues  are 
discussed.  The  simulation  consists  of  computing  the  RCS  (or 
scattered  fields)  of  a  rectangular  sheet  with  known 
resistivity.  Then  the  computed  RCS  is  used  in  the  synthesis 
equations  to  obtain  the  resistivity  function.  In  principal, 
the  synthesized  distribution  should  be  exactly  the  same  as  the 
original  distribution  used  to  compute  the  RCS.  However,  in 
practice,  the  agreement  between  the  two  depends  on  several 
factors,  such  as  the  sampling  method  used  to  fill  the  vectors 
E/  and  in  Equations  (95)  through  (98)  . 

From  the  field  equations  in  matrix  form,  the  scattered 
field  components  can  be  computed  in  any  direction  ( 6 ,  #)  when 
the  resistivity  function  Rs  is  defined.  To  form  the  vectors  Ees 
and  E^,  MxN  directions  must  be  chosen  to  compute  specific 
values.  Once  the  vectors  Ef  and  E^are  defined,  it  is  possible 
to  solve  these  equations  for  the  vectors  X a,  X*,  and  X<  The 
final  step  in  the  synthesis  procedure  is  the  reconstruction  of 
the  expansion  coefficients  {rmn}  from  these  vectors. 

One  approach  to  choosing  the  observation  directions  is  to 
define  MxN  points  in  direction  cosine  space  as  described  in 
Chapter  III-C  and  Figure  5.  The  square  region  defined  by 
-1  zu,  vzl  is  divided  into  patches  and  their  centers  represent 
the  chosen  directions.  Equally  spaced  points  in  u  and  v  have 
been  used  in  most  calculations,  but  the  spacings  in  these  two 
principal  axes  can  be  different  in  general.  The  advantage  of 
this  approach  is  that  equal  sample  spacings  work  well  with 
scattering  patterns  that  are  periodic  in  u  and  v. 

At  this  point,  two  variations  have  been  investigated. 
They  are  referred  to  as  method  (I)  and  method  (II) .  In  method 
(I)  the  observation  points  are  restricted  to  be  inside  of  the 


29 


unit  circle  (visible  region),  which  makes  physical  sense.  A 
disadvantage  of  this  method  is  the  necessity  of  forcing  the 
number  of  patches  on  the  sheet  to  be  equal  to  the  number  of 
sample  points  in  the  visible  part  of  direction  cosine  space. 
This  condition  becomes  very  restrictive  as  the  spacing  in  u 
and  v  gets  close.  Consequently,  if  MxN  points  are  to  be 
contained  in  the  unit  circle,  the  -lsu,  vsl  region  has  to  be 
divided  into  more  than  MxN  patches,  say  Nu  by  Nv.  However, 
despite  this  restriction,  the  method  was  successful,  even  for 
large  sizes  of  the  sheet. 

A  modified  version  of  this  approach  is  method  (II)  in 
which  the  MxN  direction  cosine  samples  are  spread  over  the 
complete  range  -l£U,vsl.  For  this  case  NU=N  and  NV=M.  Note 
that  the  points  outside  the  unit  circle  are  not  in  the  visible 
region  and  therefore  one  might  expect  that  these  values  would 
not  affect  the  synthesis  result.  However,  an  attempt  to  set 
the  exterior  points  equal  to  zero  results  in  a  badly  scaled 
matrix  for  the  fields,  and  the  matrix  is  almost  always 
impossible  to  invert.  The  reason  is  that  setting  the  exterior 
field  values  to  zero  is  an  attempt  to  synthesize  a  resistivity 
distribution  for  a  discontinuous  pattern.  Realizable 
resistivity  functions  do  not  yield  scattering  patterns  that 
abruptly  go  to  zero.  To  avoid  this  problem  the  calculation  can 
be  extended  to  the  exterior  points  using  complex  angles.  In 
this  case  the  expression 

cos  (0J  =  Jl ~{uk2  +  vk2)  (100) 


can  take  on  complex  values.  At  the  corners  and, 

|COS  (0*.)|  =  1 


arg{c os  (0*)  )  =  ±  . 


(101) 


While  the  first  method  makes  more  physical  sense,  the 
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second  is  more  flexible  in  choosing  the  number  of  far- field 
points  (or  patches) ,  and  is  shown  to  work  well  over  a  large 
range  of  values  of  M  and  N. 

To  verify  the  synthesis  procedure,  the  scattered  field  is 
calculated  using  two  different  resistivity  profiles:  (i)  a 
constant  resistivity  (J?s=377fi)  and  (ii)  a  two-dimensional 
linear  taper  given  by  the  equation 

Rs(x,y)  =15  0(|x|+|y|)  ft  .  (102) 

For  the  field  calculation,  a  large  enough  number  of  patches  is 
used  to  sufficiently  approximate  the  linear  shape  of  Equation 
(102)  .  The  synthesis  was  tested  using  a  few  specific  field 
points  in  the  range  of  9  to  256.  This  corresponds  to  M  and  N 
in  the  range  of  3  to  16.  Beyond  that,  the  method  still  works 
well,  but  computational  time  increases  significantly.  Roundoff 
errors  and  numerical  instability  also  affect  the  solution. 
Matrix  B  in  Equations  (92)  or  (93)  has  to  be  inverted  and  is, 
in  general,  an  ill-conditioned  matrix.  It  becomes  singular  as 
the  number  of  patches  increases.  In  such  cases  the  use  of  the 
Moore -Penrose  pseudoinverse  matrix  [Ref.  7]  rather  than  the 
inverse,  was  found  to  give  satisfactory  results  for  larger 
numbers  of  patches. 


B.  BISTATIC  CASE 

Methods  (I)  and  (II)  are  applied  for  the  bistatic  case, 
first  for  a  pure  TMj  incident  wave,  and  then  for  arbitrary 
incident  wave  polarization,  which  is  the  most  general 
situation.  For  each  resistivity  function,  the  scattered  fields 
are  computed  for  sheets  of  sizes  lXxlX  and  5Xx5X.  Data  are 
presented  for  various  numbers  of  patches  and  the  two 
resistivity  profiles.  The  plot  of  the  expansion  coefficients 
{r^}  (i.e.,  the  resistivity  profile  used  to  calculate  the 
scattered  fields)  is  compared  with  the  reconstructed 
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resistivity  function.  For  small  number  of  patches,  the 
reconstruction  is  exact  (one  can  not  distinguish  them  when  the 
plots  are  overlaid) .  But  as  the  number  of  patches  increases 
some  numerical  and  sampling  errors  can  be  observed.  In  general 
both  methods  (I)  and  (II)  work  well  in  the  entire  range  of 
synthesis  parameters  investigated. 

1.  Parallel  Polarization 

Assume  the  incident  wave  is  TM,  polarized  (  =  0  )  and  that 

it  arrives  incident  from  the  direction  6i  =  3  0°  and  4^  =  0°.  As 

an  example,  the  RCS  pattern  for  a  one  wavelength  square  sheet 
is  shown  in  Figure  6.  The  field  has  been  computed  using  method 
(II)  and  NU=NV=64 .  Various  synthesis  results  are  presented  in 
Figures  7  through  17.  Note  that 

1*  M  and  N  are  the  number  of  subdivisions  on  the 
resistive  sheet  in  the  directions  y  and  x  (see  Figure  4)  . 

2.  Nu  and  Nv  are  the  number  of  subdivisions  of  the  region 
-lsu,vsl  (see  Figure  5). 

3.  Lx  and  Ly  are  the  dimensions  of  the  sheet  in  x  and  y 
directions  expressed  in  wavelengths. 

It  should  be  noted  that  in  method  (I)  NuxNvaMxN,  while  in 
method  (II)  NU=N  and  NV=M.  In  most  cases  the  two  methods  gave 
identical  results,  thus  the  data  will  only  be  presented  once. 

In  Figures  7  and  8,  only  16  field  points  are  used  with 
constant  and  linear  taper  resistivities.  The  initial  and 
synthesized  distributions  are  identical.  Figure  9  and  10  show 
that  the  synthesized  linear  taper  resistivity  has  a  maximum 
error  of  1%  when  256  points  are  used.  The  relative  error  for 
the  constant  resistivity  case  with  the  same  number  of  points 
becomes  2%  as  shown  in  Figure  11. 

The  RCS  for  the  5\x5X  sheet  is  shown  in  Figure  12.  Again 
the  two  resistivity  functions  are  considered.  In  Figures  13 
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(constant  resistivity)  and  14  (linear  taper) ,  only  16  field 
points  are  needed  to  recover  the  exact  resistivity  function. 
In  Figures  15  through  17  it  is  seen  that  the  error  in  the 
synthesized  resistivity  is  also  negligible  when  256  points  are 
used.  It  was  found  that  for  a  sheet  with  dimensions  kXxkX,  the 
choice  of  Nu=Nv=k  gives  converged  results. 
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Figure  6:  Bistatic  RCS  of  a  resistive  sheet  with 
Lx=Ly=laxnbda,  NU=NV=64,  6=  3  0  deg,  $;=0  deg,  Rs  is  given  by 
Equation  (102) . 
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Rs  and  Ros  in  Ohms 


y  (wavelengths) 


-0.5 


x  (wavelengths) 


-0.5 


Figure  8:  Reconstructed  and  original  resistivity  functions 
are  identical.  For  this  case  Lx=Ly=lambda,  M=N=4,  N  =N=4 . 


y  (wavelengths) 


-u.o 


-u.o 


x  (wavelengths) 


Figure  9 :  Synthesized  and  original  resistivity 
Lx=Ly=lambda,  N=M=16,  NU=NV=16.  The  maximum  error  is  1 


Rs-Ros  in  Ohms 


Figure  10:  Difference  of  the  synthesized  and  original 
resistivities  shown  in  Figure  9. 
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y  (wavelengths) 


-U.3 


x  (wavelengths) 


Figure  11:  Synthesized  and  original  resistivities.  For 
this  case,  Lx=Ly=lambda,  M=N=16  and  N=NV=16.  The  error  is 


wsaP  Uj  2vepquJB|/BLu6js 


Figure  12:  Scattered  RCS  from  a  resistive  sheet.  Rs  is 
given  by  Equation  (102)  .  For  this  case  Lx=Ly=5*lambda, 
M=N=64,  NU=NV=64,  3 0  deg,  (f>=0  deg. 
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Figure  13 : 
identical . 


Synthesized  and  original  resistivity 
For  this  case,  Lx=Ly= 5* lambda,  M=N=4, 


are 
Nu=4 , 


y  (wavelengths) 


-o 


x  (wavelengths) 


Figure  14:  Synthesized  and  original  resistivity  a] 
identical.  For  this  case,  Lx=Ly=5* lambda,  M=N=4 ,  NU=NV=4 


Figure  15:  Synthesized  and  original  resistivity.  For  this 
case,  Lx=Ly=  5* lambda,  M=N=16,  NU=NV=16. 
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Figure  16:  Difference  of  the  reconstructed  and  the 
original  resistivity  distribution.  The  relative  error  is 
negligible . 
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y  (wavelengths) 


-3 


x  (wavelengths) 


Figure  17:  Original  and  synthesized  resistivity.  For  this 
case,  Lx=Ly=5* lambda,  M=N=16,  NU=NV=18.  The  error  is 
negligible . 


2.  Arbitrary  Polarization 

For  the  general  case  where  the  incident  field  has  both 
parallel  and  perpendicular  components,  both  the  E/  and  the  E/ 
vectors  given  by  Equations  (95)  and  (96)  depend  on  X,  and  X^. 
The  vectors  X^  and  Zb  are  not  independent  and  both  implicitly 
contain  the  expansion  coefficients  {rmn} .  An  expression  of 
these  equations  in  such  a  form  that  permits  the  solution  for 
(r™}  once  the  field  vectors  are  known  is  desired.  One  approach 
is  to  specify  the  magnitude  of  the  total  scattered  field  as 
given  by 


but  this  leads  to  very  complicated  functions  of  a  and  b. 
Alternately,  Equations  (95)  and  (96)  can  be  expressed  as 


<m>  *^2  • 


(103) 


Solving  this  system  for  X,  and  Xh  yields 

and 


(104) 


[-  0<t>  *-  4>0  4*4,  *-  oe]  ± 


(105) 


Once  X,  or  Xb  is  known,  Equation  (92)  or  (93)  can  be  used  to 
solve  for  the  vectors  a  or  b.  The  expansion  coefficients  can 
easily  be  determined  using  Equation  (53)  or  (65) . 

This  procedure  was  applied  to  the  resistivity  functions 
of  the  previous  section  for  sheets  of  dimensions  lXxlX  and 
5Xx5X.  Method  (II)  was  used  for  convenience  and  the  results 
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are,  for  all  practical  purposes  identical  with  the  original 
resistivity.  Figures  18  and  19  compare  the  original  and 
synthesized  resistivity  profiles  (linear  taper)  for  the  lXxlX 
sheet  and  Figures  20  and  21  for  the  5Xx5X  and  constant 
resistivity.  In  both  cases  144  points  have  been  used.  For  all 
the  cases  it  was  assumed  that  Ee'  =  E^'  =1. 
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Figure  18:  Synthesized  and  original  resistivity 
distribution  for  arbitrary  polarization.  For  this  case, 
Lx=Ly= lambda,  M=N=12,  NU=NV=12 . 
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Rs-Ros  in  Ohms 


Figure  19:  Difference  in  the  reconstructed  and  the 
original  distributions.  The  relative  error  is  less  than 
1%. 
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Rs  and  Ros  in  Ohms 


Figure  20:  Original  and  synthesized  resistivities  for 
arbitrary  polarization.  For  this  case  Lx=Ly=5*lambda, 
M=N=12 ,  NU=NV=12 ,  Es  =  =1. 
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Figure  21:  Difference  in  the  reconstructed  and  the 
original  resistivity  distributions  shown  in  Figure  20. 
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C .  MONOSTATIC  CASE 


The  synthesis  procedure  can  be  applied  to  the  monostatic 
case  using  Equation  (97)  or  (98)  for  the  scattered  field. 
Alternatively,  Equation  (86)  of  the  RCS  can  be  used,  which  is 
the  same  for  both  polarizations. 

Methods  (I)  and  (II)  have  been  investigated  for  two 
resistive  sheets  of  sizes  lXxlX  and  3Xx3X,  and  the  two 
resistivity  profiles  considered  previously;  namely  a  constant 
of  377Q,  and  a  linear  taper  given  by  Equation  (102) .  The  two 
methods  give  identical  results. 

Figure  22  compares  the  corresponding  data  for  the  lXxlX 
sheet  with  constant  resistivity.  In  Figure  23  the  difference 
between  the  synthesized  and  original  distributions  is  plotted 
for  a  linear  taper  profile  on  a  4Xx4X  sheet. 
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Rs  and  Ros  in  Ohms 


Figure  22:  Original  and  synthesized  resistivities  for  the 
monostatic  case.  On  this  case,  Lx=Ly= lambda,  M=N=12, 
N„=NV=12 . 


Rs-Ros  in  Ohms 


Figure  23:  Difference  in  the  synthesized  and  original 
linear  taper  resistivity  for  the  monostatic  case.  For  this 
case,  Lx=Ly=3* lambda,  M=N=12,  NU=NV=12. 
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D.  RELATION  TO  THE  SAMPLING  THEOREM 

In  the  previous  two  sections  it  was  demonstrated  that  the 
synthesis  procedure  works  well  for  both  methods  (I)  and  (II) 
over  a  wide  range  of  synthesis  parameters.  These  calculations 
have  verified  the  synthesis  equations  and  the  computational 
approach  to  solving  them.  So  far  the  question  of  how  many 
field  points  are  needed  to  uniquely  determine  an  unknown 
surface  resistivity  function  has  not  been  addressed. 

The  most  critical  part  of  the  synthesis  procedure  is  the 
inversion  of  the  (MxN)  by  (MxN)  matrix  B.  It  is  desirable  to 
use  the  smallest  number  (MxN)  for  computational  and  accuracy 
reasons.  To  find  the  minimum  (MxN),  a  technique  is  borrowed 
from  signal  analysis.  The  scattered  field  is  computed  using  a 
sufficiently  large  number  of  patches.  The  FFT  is  then 
performed  on  this  array  of  field  values.  The  result  can  be 
translated  into  the  minimum  number  of  patches  required  to 
accurately  represent  the  field.  This  is  based  on  the  fact  that 
plotting  a  field  pattern  (function)  is  basically  sampling  its 
distribution  adequately.  The  interval  between  the  sampled 
field  points  at  which  the  field  is  computed  is  the  sampling 
interval  and  the  reciprocal  of  this  is  the  sampling  frequency. 
Dividing  the  u  and  v  region  into  Nu  and  Nv  segments 
respectively,  the  corresponding  sampling  intervals  are  2/Nu 
and  2/Nv. 

Nyquist ' s  theorem  [Ref.  6]  requires  that  at  least  twice 
the  maximum  'frequency  component'  should  be  used  as  the 
smallest  sampling  rate.  Utilizing  the  fact  that  the  FFT  is 
symmetric  with  respect  to  the  half  of  the  sampling  frequency 
used,  i.e.,  symmetric  to  the  Nu/4  and  Nv/4  point,  a  plot  of  the 
two-dimensional  FFT  vs  the  number  of  segments  Nu  and  Nv  can  be 
obtained.  For  simplicity  an  equal  number  of  points  have  been 
used  (NU=NV)  and  the  FFT  is  performed  on  the  rows  and  the 
columns  of  the  matrix  of  the  field  values  to  identify  the 


55 


spatial  frequencies  in  u  and  v.  Figure  24  shows  the  FFT 
performed  on  the  columns  for  the  field  depicted  in  Figure  12 
using  methods  (I)  and  (II) .  As  expected,  the  frequency  content 
is  the  same  in  both  cases. 

In  Figure  25  the  spatial  frequency  content  of  the  field 
in  Figure  6  is  shown.  One  can  conclude  that  Nu  and  Nv  greater 
than  4  is  sufficient.  In  Figure  26  the  field  in  the  0=0  plane 
is  shown.  The  dotted  line  represents  the  field  using  only  16 
points  (M=N=4 )  and  the  solid  line  the  field  when  M=N=64 . 

However  for  the  5\x5\  sheet  with  the  field  shown  in 
Figure  12  and  the  frequency  content  shown  in  Figure  27,  6 
segments  per  direction  should  be  sufficient.  In  Figures  28  and 
29,  the  corresponding  field  plots  in  the  0=0  plane  show  the 
convergence . 

Figure  30  shows  the  monostatic  RCS  of  a  XxX  resistive 
sheet.  The  spatial  frequency  content  of  the  field  is  shown  in 
Figure  31,  and  it  reveals  that  4  segments  are  sufficient.  The 
corresponding  field  plots  in  Figure  32  and  33  illustrate  the 
convergence.  The  same  process  is  applied  to  the  3Xx3X  plate. 
The  frequency  content  is  shown  in  Figure  34  and  the 
corresponding  convergence  is  shown  in  Figures  35  and  36. 


56 


abs(FFT)  abs(FFT) 


Figure  24:  Frequency  content  in  u  direction  of  the 
bistatic  scattered  field  shown  in  Figure  12  using  methods 
(I)  and  (II) . 
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Sigma/lambdaA2  in  dBSm 


Angle  theta  in  deg.  (phi=0) 


Figure  26:  Bistatic  scattered  RCS  in  the  0=0  plane 
corresponds  to  Figure  6.  The  dotted  line  represent  the  RCS 
using  the  synthesized  resistivity  and  M=N=4 . 


abs(FFT)  abs(FFT) 


Angle  theta  in  deg.  (phi=0) 


Figure  29:  A  cut  on  the  u  axis  corresponding  to  Figure  12. 
Comparison  of  original  and  synthesized  patterns  for 
M=N=12 . 
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Sigma/lambdaA2  in  dBSm 
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v=sin(theta)*sin(phi) 


-1  -1 


u=sin(theta)*cos(phi) 


Figure  30:  Monostatic  RCS  from  a  one  wavelength  square 
resistive  sheet.  Rs  is  given  by  Equation  (102) . 
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Figure  33:  Using  M=N=4,  the  monostatic  RCS  due  to  the 
reconstructed  resistivity  profile  is  converged. 
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Figure  34:  FFT  of  the  monostatic  scattered  field  from  a 
resistive  sheet  with  Lx=Ly=3* lambda. 
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Figure  35:  The  monostatic  scattered  field  due  to  the 
reconstructed  resistivity  (dotted)  compared  to  the 
original  for  M=N=4. 
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E.  SUMMARY 


In  this  chapter  the  synthesis  procedure  was  computer 
simulated  and  the  corresponding  data  were  presented.  Two 
methods  were  investigated  for  choosing  observation  points  for 
generating  synthesized  distributions.  In  method  (I)  points 
inside  of  the  visible  region  were  used.  This  gives  good 
results,  but  it  was  shown  to  be  very  restrictive  because  the 
number  of  observation  points  in  the  visible  region  must  match 
the  number  of  subdomains  on  the  plate.  By  allowing  complex 
angles  it  is  possible  to  use  observation  points  in  the  entire 
region  bounded  by  -lsu,vsl  (method  (II)).  This  is  much  more 
flexible  and  still  gives  good  results. 

Using  Fourier  analysis  and  the  sampling  theorem  it  was 
possible  to  identify  the  minimum  number  of  patches  that  can  be 
used  in  the  synthesis  procedure.  In  general  the  larger  the 
sheet  the  more  points  are  necessary,  and  consequently  the 
computing  time  increases  significantly. 

To  illustrate  the  procedure  two  resistive  sheets  were 
considered,  lAxiX  and  5\x5\  for  the  bistatic  case  and  lAxiA 
and  3  Xx3  A  for  the  monostatic  case.  The  method  was  tested  first 
from  the  computational  point  of  view.  For  the  resistivity 
distributions  considered  (constant  and  linear  taper)  it  was 
found  that  patch  dimensions  of  approximately  0.5X  and  0.25X 
were  sufficient  in  the  bistatic  and  monostatic  cases, 
respectively,  to  obtain  good  agreement  with  the  original 
scattered  field. 
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V.  CONCLUSIONS 


A  method  of  RCS  synthesis  for  planar  resistive  sheets  has 
been  presented.  The  synthesis  equations  were  developed 
assuming  that  the  resistive  sheet  boundary  conditions  hold  at 
all  points  on  the  surface.  Therefore  the  equations  implicitly 
neglect  the  effect  of  the  edges  on  the  current. 

The  equations  were  verified  by  first  calculating  the  RCS 
of  the  resistive  sheet  for  a  known  resistivity  function  and 
then  using  these  values  for  back  solution  to  obtain  a 
synthesized  surface  resistivity.  In  principle  the  two 
resistivity  functions  should  be  identical,  but  because  of 
numerical  roundoff  the  agreement  is  not  exact.  However,  the 
results  for  a  wide  range  of  synthesis  parameters  were  in  very 
close  agreement  with  the  original  resistivity  distributions. 

It  was  also  shown  that  a  Fourier  decomposition  of  the 
scattering  pattern  yields  guidelines  for  choosing  the 
subdomain  size,  or  equivalently,  the  number  of  far  field 
observation  points.  It  was  found  that  a  subdomain  size  of 
about  0 . 5 A  for  bistatic  RCS  and  0.25X  for  monostatic  RCS  give 
converged  results.  By  convergence  is  meant  that  RCS  for  the 
synthesized  distribution  is  essentially  the  same  as  the  RCS 
for  the  original  scattered  distribution. 

It  is  important  to  note  that  the  synthesis  is  based  on  a 
particular  assumption  with  regard  to  the  form  of  the  current. 
The  resistive  sheet  boundary  condition  is  equivalent  to  the 
physical  optics  approximation  in  the  sense  that  the  current  in 
the  vicinity  of  the  edge  is  not  disturbed  by  the  presence  of 
the  edge.  The  current  is  only  a  function  of  the  incident  field 
and  the  surface  resistivity  (see  Equation  (4))  .  Thus  using  the 
synthesized  resistivity  distribution  in  a  method  of  moments 
code  could  yield  a  RCS  substantially  different  from  that  based 
on  the  assumed  current  if  surface  waves  exist.  Examples  are 
shown  in  Figures  37  and  38.  In  these  cases  the  surface  waves 
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are  not  as  strong  as  they  would  be  for  a  surface  with  more 
highly  conducting  edges,  and  therefore  the  agreement  between 
the  two  curves  is  relatively  close.  This  problem  is  overcome 
if  the  current  is  treated  as  an  unknown  along  with  the 
resistivity.  This  approach  allows  a  rigorous  solution  of  the 
problem  and  is  recommended  that  future  work  is  directed  along 
these  lines. 
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Figure  37:  Comparison  of  the  monostatic  RCS  using  MM  and 
the  approximate  method.  For  this  case,  Lx=Ly= lambda, 
M=N=12 ,  NU=NV=12  .  R,  =  377Q. 
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Figure  38:  Comparison  of  the  monostatic  RCS  using  MM  and 
the  approximate  method.  Lx=Ly=lambda ,  M=N=12,  NU=NV=12.  Rs 
is  given  by  Equation  (102)  . 
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APPENDIX.  COMPUTER  CODES 


This  appendix  contains  the  codes  developed  to  computer 
simulate  the  synthesis  equations.  Code  field. m  calculates  the 
scattered  field  or  RCS  for  each  case  (monostatic,  bistatic, 
method  (I),  method  (II)),  given  the  size  of  the  resistive 
sheet  and  the  surface  resistivity.  Code  synth. m  uses  the  field 
values  calculated  from  field. m,  to  synthesize  the  surface 
resistivity.  Code  synthld.m  uses  the  surface  resistivity 
values  calculated  by  synth. m  to  plot  the  RCS  in  the  0=0  plane. 
Code  fanalysl.m  performs  the  FFT  analysis  for  methods  (I)  and 
(II),  and  code  fanalys2.m  performs  the  FFT  analysis  in  the  u 
and  v  direction  of  the  scattered  field.  Function  coord. m  is 
used  in  the  first  two  codes  to  find  the  coordinated  of  the 
patches  given  the  size  of  the  sheet  and  the  number  of  patches. 
Finally,  program  fit.m  finds  the  closest  number  of  points 
closest  to  an  input  value  that  fit  in  the  unit  circle,  and  is 
used  in  method  (I) . 
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Program  field. m 


%  This  program  calculates  for  a  rectangular  sheet: 

%  (i)  Bistatic  case:  (sigma_theta_theta,  or 
sigma_phi_phi) /lambdaA2 


%  or  the 

%  scattered  field  for  the  general  case 

of  the 


Eth  i  and 


incident  field,  given  the  values 


%  Ephi_i  and  assuming  Co=l. 

%  (ii)  Monostatic  case:  sigma/lambdaA2 


%  The  plate  can  be  PEC  or  resistive.  Plots  of  the 

resistance, 

%  as  well  as,  3-D  plots  of  the  field  are  shown. 

%  Works  in  the  direction  cosine  space. 


%  Date:  10  Nov  1994 
%  By:  Nick  Faros. 


clear 


S' S' S' S' S' S- 
15  o  o 


INPUT  PARAMETERS 


o  o  o  o  o  o 
o  o  o  o  o  o 


choicel=input (' Enter  1  for  bistatic  or  0  for  monostatic:  '); 
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choice2=input ( 'Enter  1  for  method  (II)  or  0  for  method  (I): 

'  )  ; 

al=input (' Enter  size  in  -x  (multiples  of  wavelength):  '); 
bl=input (' Enter  size  in  -y  (multiples  of  wavelength):  '); 

%  Assume  wavelength  1=1 
1  =  1; 

k=2*pi/l;  %  wavenumber 

a=al*l;  %  dimension  in  x  (in  meters) 

b=bl*l;  %  dimension  in  y  (in  meters) 

Nx= input (' Enter  #  of  patches  in  x  (N) :  '); 

Ny=input (' Enter  #  of  patches  in  y  (M) :  '); 

dx=a/Nx;  dy=b/Ny; 

xx=-a/2 : a/ (Nx- 1) :a/2;  yy=-b/2 :b/ (Ny- 1) :b/2; 

%  xx, yy  are  the  dimension  of  the  plate  in  meters 


R=input ( 'Resistance: 0=PEC,  l=constant,  2=linear  taper, 
3=ramp : ' ) ; 
if  R==l 

cl=input (' Enter  constant  value:  '); 

end 
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xl=reshape (xo, 1 , Nx*Ny) ; 
yl=reshape (yo, l,Nx*Ny) ; 


o  o  o  o  o  o 


FIND  PATTERN  OF  THE  INPUT  RESISTANCE 


Some  cases 


if  R=  =  0 

r=zeros (Ny , Nx) ;  %  Zero  everywhere 

elseif  R==l 

r=cl*ones (Ny , Nx) ;  %  Constant 

elseif  R==2 

r=abs (xo) +abs (yo) ;  %  linear  taper  in  both  dimensions 

r=150*r ; 
elseif  R==3 

r=zeros ( Ny , Nx ) ; 
for  1=1 : Nx 

r  ( :  , 1) =100*l*ones (Ny, 1) ;  %  ramp 

end 

end 

rl=reshape (r, 1 , Nx*Ny) ; 


O  O  O  O  O  O  O  O  O  O  oo^^^^^'^'^'o^’o^'o'o'o'o'o'o’o^'o'o'o'o'o'o 


5'5'5'S-5-5-5-5'5'5-$-9-$-5--S-9-0-5-o,o,o^o,o,o^o, 


%  Next  find  coordinates  of  the  center  of  the  patches  in 
%  direction  cosine  space.  If  method  (I)  is  used,  run  fit.m 
%  first,  to  find  nx,  ny.  If  method  (II)  is  used  nx=Nx  and 
%  ny=Ny; 


if  choice2==l  %  Method  (II) 
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nx=Nx; 
ny=Ny ; 


elseif  choice2==0  %  Method  (I) 
nx=4  %  Nu 

ny=5  %  Nv 

end 


[u,v] =coord (2,2, nx, ny) ; 
ul=reshape (u, 1 , nx*ny) ; 
vl=reshape (v, 1 , nx*ny) ; 


S-S'S-S'S'S'S'S'S'S-S'S'S'S'S'S'S'S^S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S^S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S'S' 

'o'o’o^'o'o^'o^'o^'o^'o'o'o'o'o'o'o'o'o'o'o^'o^'o^'o^o'o'o^^^^^^o^oo^^^o^o^^^^^oo^^^o^ooo^^o 


sintheta=sqrt (ul . A2  +vl.A2); 

I=f ind (sintheta>l) ; 
phi=atan2 (vl,ul) ; 
costheta=sqrt (1-sintheta. A2) ; 
if  choice2==0  %  Method  (I) 

phi(I)  =  [];  costheta  (I)  =  []  ;  sintheta  (I)  =  []  ; 

ul  (I)  =  []  ;  vl(I)  =  []  ; 

end 


2- S' S' S' S' S' S' S' S' S' S' S' S' S' S' S' S' S' 2- S' S' S' S' S' S' S' S' S' S' 2- 2- 


,  q.  p,  q,  Q,  o,  q. 


'Q-Q,q,q,q,q,q,q,q.q,g,q.q,o^q, 


"o  "o  0'b‘o‘b'b  o'b'o'b  o  oo’b'ooooo  o  oooo 


O  O  O  O  O  0-0-0- 

oooooooo 


BISTATIC  CASE 


S' S' S' S' S' S' S' S' S' S' 


if  choicel==l  %  Bistatic 

thid=input (' Enter  theta  inc.  in  degrees:  '); 
thir=thid*pi/180 ;  %  translate  it  into  radians 

phiid=input (' Enter  phi  inc.  in  degrees:  '); 
phiir=phiid*pi/180 ; 
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coef 1  =  1 . / (2*rl+377*cos (thir)  )  ; 
coef 2  =  1 . / (2*rl*cos (thir) +377)  ; 


bistatic 

bistatic 


o, 
o 

% 

%  u=sin (theta) *cos (phi) 

%  v=sin (theta) *sin (phi) 

ui=sin (thir) *cos (phiir) ; 
vi  =  sin (thir) *sin (phiir)  ; 


sincl=sin (k/2*dx* (u (1, : ) +ui) ) ./ (k/2*dx* (u (1, : ) +ui) ) ; 
Sl=f ind (isnan (sincl) ) ; 
sincl (SI) =ones (size (SI) ) ; 

sinc2=sin(k/2*dy* (v ( : , 1) +vi) ) ./ (k/2*dy* (v( : , 1) +vi) ) ; 
S2=f ind (isnan (sinc2 ) ) ; 
sinc2 (S2 ) =ones (size  (S2 ) ) ; 

sinc=sinc2* sincl ; 
sinc=reshape (sine, l,nx*ny) ; 
if  choice2==0 
sine  (I)  =  []  ; 
end 


choice3=input  (' Enter  1  for  sigxna_theta_theta ,  2  for 
sigma _phi_phi,  or  3  for  the  general  scattered  field:  '); 

if  choice3==l  %  Compute  sigma_theta_theta 

for  p=l : Ny*Nx 

expn=exp ( j*k* (xl* (ul (p) +ui) +yl* (vl (p) +vi) ) ) • 
si (p)  =  (expn*conj (coef 1 ' )  )  ; 
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end 


%E=dx*dy*cos ( thir) * cos theta . *sinc . *cos (phi-phiir) . *sl ; 
E=dx*dy*cos (thir) *costheta . *sinc . *cos (phiir) . *sl ; 


RCS1= ( (377*2*pi*E) . x2) ./pi; 

RCS=10*logl0 (abs (RCS1) ) ; 

RCS=reshape (RCS,Ny,Nx) ; 

elseif  choice3==2  %  Compute  sigma_phi_phi 

for  p=l:Ny*Nx 

expn=exp ( j  *k* (xl* (ul (p) +ui) +yl* (vl (p) +vi) ) ) ; 
s2 (p)  =  (expn*conj (coef2 ' )  )  ; 
end 

E=dx*dy*cos (thir) *sinc . *cos (phi-phiir) . *s2 ; 

RCS1= (377*2*pi*E) .^2. /pi; 

RCS=10*logl0 (abs (RCS1) ) ; 

RCS=reshape (RCS,Ny,Nx) ; 

elseif  choice3==3  %  Compute  field 


Eth_i  =input (' Enter  magnitude  Eth_i:  '); 
Ephi_i  =input (' Enter  magnitude  Ephi_i :  '); 

for  p=l:Ny*Nx 

expn=exp ( j  *k* (xl* (ul (p) +ui) +yl* (vl (p) +vi) ) ) ; 
si (p)  =  (expn*conj (coef 1 '  )  )  ; 
end 
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for  p=l:Ny*Nx 

expn=exp ( j  *k* (xl* (ul (p) +ui) +yl* (vl (p) +vi) ) ) ; 
s2 (p) = (expn*conj (coef 2 ' ) ) ; 
end 


E_th=dx*dy*cos ( thir) * cos theta . *sinc . * (Eth_i*cos (phi -phi ir) . * 
sl+Ephi_i*sin (phi -phiir )  .*s2}  ; 


E_phi=-dx*dy*cos (thir) *sinc. * (Eth_i*sin (phi-phiir) . *sl-Ephi 
i*cos (phi-phiir) .  *s2 ) ; 

E=sqrt ( (abs (E_th) ) .^2+ (abs (E_phi) ) .^2) ; 

El=reshape (E, Ny , Nx) ; 

end 

if  choice3==3 

save  fnb  a  b  E_th  E_phi  Eth_i  Ephi_i  r  ui  vi  thir 
phiir  ul  vl  xl  yl  xx  yy  choice3  choice2  %  Use  in  synth. m 
else 

save  fnb  a  b  E  r  ui  vi  thir  phiir  ul  vl  xl  yl  xx  yy 
choice3  choice2  %  Use  in  synth. m 
end 


oooooooooo 


MONOSTATIC  CASE 


%% 


S-S-S-S-5-S'S'S-S- 


elseif  choicel==0 

sincl=sin (k*dx*u {!,:))./ (k*dx*u (1, : ) ) ; 
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Sl  =  f ind (isnan (sincl)  )  ; 
sincl (SI) =ones (size (SI) ) ; 

sinc2=sin (k*dy*v( : , 1) ) . / (k*dy*v ( : , 1)  )  ; 
S2=find( isnan (sinc2) ) ; 
sinc2 (S2 ) =ones (size (S2)); 

sinc=sinc2* sincl; 
sinc=reshape (sine, l,nx*ny) ; 

coef=l . / (2*rl+377) ;  %  monostatic 

for  p=l:Ny*Nx 

expn=exp (2* j *k* (xl*ul (p) +yl*vl (p) ) ) ; 
s (p) = (expn*conj (coef ' ) ) ; 
end 

E=dx*dy*sinc . *s ; 

RCS1= (377*2*pi*E) .A2/pi; 

RCS=10*logl0 (abs (RCS1) ) ; 

RCS=reshape (RCS , Ny , Nx) ; 

if  choice2==l  %  Method  II 


save  fnm  a  b  E  ul  vl  xl  yl  r  xx  yy  choice2  %  Use  it  in 
synth.  m 

elseif  choice2==0  %  Method  (I) 

E  ( I )  =  [  ]  ; 

save  fnm  a  b  E  ul  vl  xl  yl  r  xx  yy  choice2  %  Use  it  in 
synth  .m 

end 

end 
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save  nf  E  Nx  Ny  I 


Use  in  fanalys.m 


time=toc 


9-  2-  2-  2-  2  2  2  2  2-  DT  DTC  9-  2*  2  2  2  5-  9-  9-  2 

'o'o'o'o'o'o'o'o'o  iri-iUlo  *b  *0  “b  *b  *o  *6  o  o  o  o 

if  choicel==l 

if  choice2==l  %  Method  (II) 
if  choice3==l  j  choice3==2 

figure (1) ;mesh (xx,yy, real (r) ), xlabel (' Size  in  x' ) 
ylabel('Size  in  y ' ) , zlabel ( ' Ohms ' ) 


figure (2 ) 
mesh (u, v, RCS) 

xlabel ( ' u=s in (theta) *cos (phi) ' ) 
ylabel ( ' v=s in (theta) *sin (phi) ' ) 
zlabel ( ' Sigma/lambda^2  in  dBSm' ) 

sint=u (1 ,  :  )  ; 

th=asin (sint ) ;  th=th*180/pi ; 
figure (3 ) 

plot ( th, RCS (Ny/2 , : ) ) , grid 

xlabel ( 'Angle  theta  in  deg.  at  phi=0  plane') 
ylabel (' Sigma/ lambda^ 2  in  dBSm') 

RCSl=reshape (RCS1 , Nx, Ny) ; 

mesh (u, v, abs (RCS1) ) 

end 

end 

elseif  choicel==0 

figure (1) ;mesh (xx,yy, real (r) ), xlabel (' Size  in  x' ) 
ylabel ('Size  in  y '), zlabel (' Ohms ' ) 
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figure (2) 
mesh (u, v, RCS) 

xlabel ('u=sin (theta) *cos (phi) ' ) 
ylabel ( 'v=sin (theta) *sin(phi) ' ) 
zlabel ( ' Sigma/lambdaA2  in  dBSm' ) 

sint=u (1, : ) ; 

th=asin (sint ) ;  th=th*180/pi ; 
figure (3) 

plot (th, RCS (Ny/2, :) ) ,grid 

xlabel ( 'Angle  theta  in  deg.  at  phi=0  plane') 
ylabel (' Sigma/lambdaA2  in  dBSm') 

RCSl=reshape (RCS1 , Nx, Ny) ; 

mesh (u, v, abs (RCS1) ) 

end 

end 

end 

save  mf  RCS  th 
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Code  synth. m 


o, 
o 

-o  This  code  reads  from  the  files  fn*.mat  (output  of 
■°  the  code  field. m)  ,  RCS  or  field  values  at  specific 
%  points,  and  synthesizes  the  surface  resistivity. 

%  This  code  divides  the  sheet  and  DCS  into  'big' 

■z  patches  compared  to  the  field. m  code.  The  number  of 
%  patches  should  be  choosen  such  that  the  center  of 
%  the  big  patch,  coincides  with  the  center  of  a  small. 

•6  So,  the  big  patch  must  contain  an  odd  number  of 
%  smaller  patches. 

%  Note  that  the  codes  field. m  and  synth. m  can  also 
%  have  the  same  number  of  patches. 

%  Works  for  monostatic  and  bistatic  case, 
clear 

choicel=input ('Enter  1  for  bistatic  or  0  for  monostatic:  '); 

if  choicel==l 
load  fnb 

elseif  choicel==0 
load  fnm 

end 

k=2*pi ; 

Nx=input (' Enter  #  of  patches  in  x  (N)  :  ')■ 

Ny=input (' Enter  #  of  patches  in  y  (M) :  '); 

%  Nx,  Ny  are  the  'big'  patches, 
tic 

Dx=a/Nx;  Dy=b/Ny; 

XX=-a/2:a/ (Nx- 1) :a/2;  YY=-b/2:b/ (Ny-1) :b/2; 
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2-2-2-2-2- 
o  o  o  ^  ^ 


FIND  COORDINATES  OF  THE  PATCHES  OF  THE  SHEET 


2-  2-  2-2-2- 
^  o  o  ^  ^ 


[xo, yo] =coord (a, b, Nx, Ny) ;  %  Call  function  coord. m 

xll=reshape (xo; l,Nx*Ny) ; 
yll=reshape (yo, 1 , Nx*Ny) ; 


SrSr§r?r?r§r?r?r%?r?r?r9'e's- 

oooooooooooo 


^^O-O.q.q.O.CVOOO.OOOOO.O.O.O.O.O.O.O.O.O^O.O.O.O.^O.q.O.O.O.O.O.O.O.O.O.O.O.CUO.O.O. 


2-2-2-%%% 
o  o  o  o  o  o 


FIND  COORDINATES  OF  THE  PATCHES  IN  DCS 


if  choice2==l 
nx=Nx ; 
ny =Ny ; 


[U, V] =coord (2 , 2 , nx, ny) ; 
Ul=reshape (U, 1 , nx*ny) ; 
Vl=reshape (V, 1, nx*ny) ; 


elseif  choice2==0 
nx=4  %  Nu 

ny=5  %  Nv 

[U, V] =coord(2,2,nx,ny) ; 
Ul=reshape (U, 1, nx*ny) ; 
Vl=reshape (V, l,nx*ny) ; 
sintheta=sqrt (U1 . a2  +V1.a2); 

II=f ind (sintheta>l)  ; 

U1(II)  =  []  ;  V1(II)  =  []  ; 
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end 


J-5'^5'S'?'S'S'?'J'2'5'9'5'&'^S'S'^5'2'S'S'5'2'S'S'S'5'S'5'2'5-^S-^5-2'S'^2-5,5<5,2,S'5-5'2'^5.5-S,5-9,9.o,q,q,o, 

o  o  o  o  o^^^^oooooo^oooo^^^^^^^^^^^'o-o^-o^'o'o'o'o'o'o'o'o-o'o-o^'o'o^'o^'o'o'o'o'o'o'o'o'o'o'o'^ 

%-  £•  ^  5-  2-  ^  £- 

T>  O  O  O  O  O  T> 


o  oo  on 
o  o  o  "o  "o 


This  part  finds  the  index  of 


the  selected  points 


&•  s- e- s- e- e- % 

o  o  o  o  o  o  o 


for  1=1 : length (Ul) 
f ll=ul -U1 (1) ; 

f 21=f ind (abs (f 11) -min (abs (f 11) ) =  =  0)  ; 
f 12=vl - VI (1) ; 

f 22  =  f ind (abs (f 12 ) -min (abs (f 12 ) ) =  =  0)  ; 


11  =  1; 

f f  =  f 22  - f 21  (11)  ; 
f f l=f ind (f f ==0) ; 
while  f  f  1  =  =  [] 
11=11+1; 
f f  =  f 22  - f 21 (11)  ; 
f f l=f ind (f f==0) ; 

end 

index (1) =f21 (11)  ; 

end 


o  o 
<r 
o  o 


o  o  o  o  o  o  o  o  o  o  o  o  o^lj^lj^j^'o^'o'o'o'o'o'o'o'o-cr-o'o'o'o'o'o'o'o'o'o'o'o'o'o-o'o 


O  O  O 


%%%%%%% 

UUl=ul (index) ; 

Wl=vl  (index)  ; 

[UU1 '  Ul '  W1 '  VI '  ]  ; 
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Xl=xl (index) ; Yl=yl (index) ; 

[XI'  xll '  Yl*  yll' ] ; 

sintheta=sqrt (U1 . a2  +V1.a2); 
pphi=atan2 (V1,U1)  ; 

ccos theta=sqrt ( 1 - sintheta . A2 ) ; 


if  choicel==l  %  Bistatic 

sincl=sin (k/2*Dx* (U(lf : ) +ui) ) ./ (k/2*Dx* (U(l, : ) +ui) ) ; 
Sl=find (isnan (sincl) ) ; 
sincl (SI) =ones (size (SI)  )  ; 

sinc2=sin(k/2*Dy* (V( : , 1) +vi) ) ./ (k/2*Dy* (V ( : , 1) +vi) ) ; 
S2=find(isnan(sinc2) ) ; 
sinc2 (S2 ) =ones (size (S2)  )  ; 

ssinc=sinc2*sincl; 
ssinc=reshape (ssinc, 1 , nx*ny) ; 
if  choice2==0 
ssinc  (II)  =  []  ; 
end 

C=  []  ; 

for  l=l:Nx*Ny 

e=exp ( j  *k* (XI* (U1 (1) +ui) +Y1* (VI (1) +vi) ) ) ; 

C=(C;e]  ; 
end 


if  choice3==l  %  Sigma_th_th 


EE=E (index) ; 
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E1=EE . / (Dx*Dy*cos ( thir) *ccostheta . *ssinc . *cos (pphi -phiir) ) ; 
E1=EE . / (Dx*Dy*cos (thir) *ccos theta . *ssinc . *cos (phiir) ) ; 


Z=f ind (isnan (El) ) ; 

El (Z) =zeros (size (Z)  )  ; 

El=conj (El ' ) ; 

%g=C\El; 
g=pinv (C) *E1 ; 

imp= . 5* (l./g-377*cos (thir) ) ; 
impl=reshape (imp,Ny,Nx) ; 

dd=real (imp) ; 

elseif  choice3==2  %sigma_phi_phi 


EE=E ( index) ; 

E1=EE . / (Dx*Dy*cos (thir) *ssinc. *cos (pphi-phiir) ) ; 
Z=f ind (isnan (El) ) ; 

El (Z) =zeros (size (Z) ) ; 

El=conj  (El ' ) ; 
g=C\El; 

%g=pinv (C) *E1; 

imp= (1/ (2*cos (thir) ) ) * (l./g-377) ; 
impl=reshape (imp,Ny,Nx) ; 

dd=real ( imp) ; 

elseif  choice3==3 
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El_th=E_th . / (Dx*Dy*cos (thir) *ccostheta . *ssinc) ; 

Zl  =  f ind (isnan (El_th) ) ; 

El_th (Zl) =zeros (size (Zl) )  ; 

El_th=conj (El_th' ) ; 

El_phi=E_phi  .  /  ( -Dx*Dy*cos  (thir)  .  *ssinc)  ; 

Z2=f ind (isnan (El_phi) ) ; 

El_phi  ( Z2 )  =zeros  (size  (Z2)  )  ; 

El_phi=conj (El_phi' ) ; 

El-  [El_th;El_phi]  ; 

all=conj (Eth_i*cos (pphi-phiir) ) ' ; 
al2=conj (Ephi_i*sin (pphi-phiir) ) ' ; 
a21=conj  (Eth_i*sin (pphi-phiir) )  '  ; 
a22=conj ( -Ephi_i*cos (pphi-phiir) ) ' ; 

Yl= ( (al2.*El_phi) - (a22.*El_th) ) ./( (al2 . *a21-all. *a22) ) ; 
%Y2= (El_th-all . *Y1) . / al2 ; 

gl=pinv (C) *Y1; 

%g2=pinv (C) *Y2 ; 


imp= . 5* (1 . /gl- 3 77* cos (thir) ) ; 
%imp=l/ (2*cos (thir) ) * (l./g2-377) • 

impl=reshape (imp,Ny,Nx) ; 
dd=real (imp) ; 


end 

save  faros  dd  Dx  Dy  k  xll  choicel  choice3  thir  phiir  ui 


'o'o'o'o'o'o'o'o'o'o'o^'o’o^'o^^o  IviUi\Uo  iili  XL  'S'b'o'o'o'b'b'b'b'b'b'S'^o'oooo 
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elseif  choicel==0 


Monostatic 


sincl=sin (k*Dx*U (1, : ) ) . / (k*Dx*U (1,  :  )  )  ; 
Sl=f ind (isnan (sincl) ) ; 
sincl (SI) =ones (size (SI) ) ; 

sinc2=sin (k*Dy*V ( : , 1) ) . / (k*Dy*V ( : , 1) ) ; 
S2=f ind (isnan (sinc2 ) ) ; 
sinc2 ( S2 ) =ones (size (S2 )  )  ; 

ssinc=sinc2* sincl  ; 
ssinc=reshape (ssinc , 1 , nx*ny ) ; 


C=  [  ]  ; 

for  l=l:Nx*Ny 

e  =  exp (2* j  *k* (X1*U1 (1) +Y1*V1 (1)  )  )  ; 
C= [ C ; e ]  ; 
end 

EE=E ( index) ; 

E1=EE . / (Dx*Dy*ssinc)  ; 

Z=f ind (isnan (El)  )  ; 

El (Z) =zeros (size (Z) ) ; 

El=conj (El ' ) ; 
g=C\El; 

%g=pinv (C) *E1; 
imp= . 5* ( 1 . /g- 377 ) ; 
impl=reshape (imp, Ny , Nx) ; 

dd=real (imp) ; 


save  faros  dd  Dx  Dy  k  xll  choicel  %  Use  in  synthld.m 
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time=toc 


end 


figure 

mesh (XX, YY, real (impl) ) 
hold  on 
mesh (xx, yy , r) 
hold  off 

xlabel ( 'x  (wavelengths) ' ) 

ylabel ( 'y  (wavelengths) ' ) 

zlabel ( ' Rs  and  Ros  in  Ohms ' ) 

rerror=max (max (abs (impl) -r) ) /max (max (r) ) 
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Code  synthld.m 


%  This  code  reads  the  surface  resistivity  from  the  file 
faros . dat 

%  created  by  field. m  and  plots  the  scattered  RCS  in  the 
principal 

%  plane  phi=0  (u  axis) . 
clear 

load  faros 

sinth=0 :  ( 1/49 )  : 1 ;  %  Sampling  the  sin  (or  the  dcs) 
costh=sqrt ( 1 - sinth .*2) ; 

if  choicel==l  %  Bistatic 

if  choice3  =  =  l  %  Sigma__th_th 

d=l . /  (2*dd+377*cos (thir)  )  ; 
for  p=l:50 

expn=exp ( j  *k* (xll* (sinth (p) +ui) ) ) ; 
s (p) = (expn*conj (d) ) ; 

end 

sincl=sin (k/2*Dx* (sinth+ui) ) . / (k/2*Dx* (sinth+ui) ) ; 
Sl=f ind ( isnan (sincl ) ) ; 
sincl (SI) =ones (size (SI) )  ; 


RCS= (377*2*pi*Dx*Dy*cos (thir) *costh*cos (phiir) . *sincl.*s) . A2 
/pi; 

elseif  choice3==2  %  Sigma  phi  phi 
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d=l . /  (2*dd*cos (thir) +377)  ; 
for  p=l:50 

expn=exp ( j  *k* (xll* (sinth (p) +ui) ) ) ; 
s (p) = (expn*conj (d) ) ; 

end 

sincl=sin (k/2*Dx* (sinth+ui) ) . / (k/2*Dx* (sinth+ui) ) ; 
Sl=f ind (isnan (sincl) ) ; 
sincl (SI) =ones (size (SI) ) ; 


RCS= (377*2*pi*Dx*Dy*cos (thir) *cos (phiir) *sincl . *s) . ^2/pi; 


elseif  choice3==3 
brake 


end 


elseif  choicel==0  %  Monostatic 

d=l . / (2*dd+377) ; 
for  p=l:50 

expn=exp (2* j  *k* (xll*sinth (p) ) ) ; 
s (p) = (expn*conj (d) ) ; 

end 


sincl=sin (k*Dx*sinth) . / (k*Dx*sinth) ; 
Sl=f ind (isnan (sincl) ) ; 
sincl (SI) =ones (size (SI) ) ; 
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RCS= (377*2*pi*Dx*Dy*sincl . *s) . ^2 /pi ; 


end 


RCSl=f liplr (RCS ) ; 

RCS= [RCS1  RCS]; 

theta=asin (sinth) ; 
thetal=f liplr (theta) ; 
theta2= [ - thetal  theta]; 
theta3=theta2*180/pi ; 
sinthetal=f liplr (sinth) ; 
sintheta2= [ - sinthetal  sinth]; 


figure ( 1 ) 

plot (theta3 , 10*logi0 (abs (RCS) ) , ' r' ) 
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%  Code  fanalysl.m 

%  This  code  calculates  the  FFT  of  the  field 
%  of  methods  (I)  and  (II) 

clear 
load  nf 
%  Method  (I) 
x=  length (E)  ; 

Yl-fft (E) ; 

MM=  [Nx  Ny]  ;MMM=min(MM)  ; 
fl=MMM/2* (0 :x/2-l) /x; 

%  Method  (II) 

E2=E ; 

E2  (I)  =  []  ; 
x2=length (E2) ; 

Y-fft  (E2 )  ; 

f 2=MMM/ 2  * ( 0 : x2 / 2  - 1 ) / x2 ; 


figure, subplot (211) 
plot (f2,abs ( Y ( 1 : x2 / 2 ) ) ) 

%axis([0  3  0  1500]); 

xlabel ( ' Nu  or  Nv' ) ;ylabel ( ' abs (FFT) ' ) 
subplot (212 ) 

plot (fl, (abs (Y1 (l:x/2) ))) ;  %axis([0  3  0  1500]); 
xlabel ('Nu  or  Nv' ) ;ylabel ( ' abs (FFT) ' ) 
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%  Code  f analys2 . , 

%  This  code  calculates  the  FFT  of  the  field 
%  in  u  and  v  directions 

clear 
load  nf 

%  Method  (I)  columns 
%E=  (377*2*pi*E)  .A2/pi; 
x=length (E)  ; 

Yl=abs (f f t (E) ) ; 

MM= [Nx  Ny] ;MMM=min(MM) ; 
fl=MMM/2* (0:x/2-l)/x; 

%  Method  (II)  -rows 
E2  =  reshape  (E,  Ny,  Ny)  ; 

E2=conj (E ' ) ; 

E2=reshape (E2, l,Nx*Ny) ; 

Y2=abs (fft ( E2 ) ) ; 

figure, subplot (211) 

plot (f 1, Y1 (1 :x/2) ) 

axis ( [0  5  0  . 8] ) ; 

xlabel ( '  Nv' ) ;ylabel ( ' abs (FFT) ' ) 

subplot (212 ) 
plot (fl, Y1 ( 1 : x/2 ) ) ; 
axis (  [0  5  0  . 8] )  ; 

xlabel ( ' Nu ' ) ; ylabel ( ' abs ( FFT) ' ) 
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function  [xo,yo] =coord(a,b,Nx,Ny) 

%  Given  the  dimensions  (a)  in  x,  (b)  in  y  of  a 
%  rectangle  and  the  number  of  segments  (Nx)  in  x 
%  and  (Ny)  in  y,  the  coordinates  of  center  of 
%  each  patch  is  returned. 

%  Point  0  of  the  coordinate  system  is  at  the 
%  center  of  the  rectangle. 

%  Date:  1  Nov  1994 
%  By:  Nick  Faros 


dx=a/Nx; 

dy=b/Ny; 

if  rem(Nx,2)==0  %  i.e.  Nx  is  an  even  number 

for  n=l:Nx 
for  i=l:Ny 

xo(i,n) =- (Nx/ 2  - (n-1) -0.5) *dx; 
end 
end 

else  %  if  Nx  is  an  odd  number 

for  n=l:Nx 
for  i=l:Ny 

xo (i , n) =- ( (Nx- 1) /2- (n-1) -0.5) *dx-dx/2 ; 
end 
end 

end 

if  rem(Ny,2)==0  %  i.e.  Ny  is  an  even  number 

for  n=l:Nx 
for  i=l:Ny 

yo (i , n) =  (Ny/2- (i-1) -0.5) *dy; 
end 
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end 


else  %  if  Nx  is  an  odd  number 

for  n=l:Nx 
for  i=l:Ny 

yo  (i  ,  n)  =  (  (Ny+1)  /2  -  (i-1)  -  0 . 5 ) *dy- dy/2 ; 
end 
end 

end 
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%  Code  fit.m 

%  This  program  fits  M  points  in  the  unit  circle. 

%  The  spacing  is  uniform  in  one  or  both  directions. 

%  Date:  1  Nov  1994 
%  By:  Nick  Faros, 
clear 

M= input (' Enter  number  of  points  in  the  unit  circle:  '); 

nx=2 ;  ny=2;  %  Start  at  those  values. 

dx=2/nx; 

dy=2/ny; 


Find  the  coordinates  of  the  center 
of  the  patches. 

[ul , vl] =coord(2,2,nx,ny) ; 
ul=reshape (ul , 1 , nx*ny ) ; 
vl=reshape (vl , 1 , nx*ny ) ; 
sintheta=sqrt (ul . a2  +v1.a2); 

I=f ind (sintheta>l) ; 
vl(I)  =  []  ;ul  (I)  =  []  ; 

L=length (vl) ; 


while  L<M 

if  rem(k,2)==0 
nx=nx+l  ; 

tul,vl] =coord(2,2,nx,ny) ; 
ul=reshape (ul , 1 , nx*ny) ; 
vl=reshape (vl , 1 , nx*ny ) ; 
sintheta=sqrt (ul . a2  +v1.a2); 
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I=f ind (sintheta>l) ; 
vl(I)  =  []  ;ul  (I)  =  []  ; 

L=length (vl) ; 

elseif  rem(k,2)~=0 
ny=ny+l ; 

tul , vl ] =coord (2 , 2 , nx,  ny) ; 
ul=reshape (ul,l, nx*ny) ; 
vl=reshape (vl , 1 , nx*ny) ; 
sintheta=sqrt (ul . A2  +vl.^2); 

I=f ind (sintheta>l)  ; 
vl(I)  =  []  ; ul  ( I )  =  []  ; 

L=length (vl) ; 

end 
k=k+l ; 
end 

[nx  ny  L] 

plot ( ul , vl , '  * ' ) , xlabel ('u' ) ,ylabel ('v' ) 
title (' Points  inside  the  unit  circle') 
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